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We generalize the usual abelian Berry phase generated for example in a system with two non- 
degenerate states to the case of a system with two doubly degenerate energy eigenspaces. The 
parametric manifold describing the space of states of the first case is formally given by the G(2, 1) 
Grassmannian manifold, while for the generalized system it is given by the G(4, 2) one. For the 
latter manifold which exhibits a much richer structure than its abelian counterpart we calculate the 
' connection components, the field strength and the associated geometrical phases that evolve non- 

I trivially both of the degenerate eigenspaces. A simple atomic model is proposed for their physical 

. implementation. 

(N ■ 

;_l I PACS numbers: 03.65.Vf, 03.67.-a, 42.50.Vk 
^ '. 

! I. INTRODUCTION 

' Geometrical phases have attracted much interest since the seminal work by Berry [1]. The simple example of 

^ ' the abeUan Berry phase produced, e.g. by the adiabatic transition of a spin-1/2 particle that follows a rotating 
magnetic field, has found many applications in quantum optics, molecular physics and so on. Theoretically, it was 

] extended to non-abelian phases by Wilczek and Zee [2] by additionally employing to the setup of the usual Berry 

^ ■ phase a degenerate structure that allows geometrical unitary evolutions of higher dimensionality describing transitions 

' of population within the degenerate eigenspace. There are different applications of geometrical evolutions in the 
literature [3] and in particular related to quantum information. Special attention has been given to the evolutions 

' that are parameterized by the n-dimensional complex projective manifold, CP". It corresponds to the parametric 

■ manifold of a physical model where n states out of n + 1 are preserved degenerate at all times [4] . For this case 

f*) i' the parametric space A4 is given by 

S ■ CP- « + ^) 

C3 ; U{n) X U{1) ' 

^3 . 

O^' where the dictated control transformations are between the non-degenerate state and each degenerate one. By 
^ , performing adiabatically closed paths in this parametric space geometrical unitary operators are generated, called 
holonomies, P, acting solely on the degenerate states. Their relevance for performing quantum computation was first 
demonstrated in [5]. Since then they enjoyed theoretical applications in quantum optical models with laser beams [6], 
trapped ions [7], optical cavities [8] or quantum dots [9]. 

A further generalization of the parametric control space is realized by employing the Grassmannian manifolds. 
They are given by the coset space structure 

G[m, n) ~ 



U{n) X U{m — n) 



The complex projective manifolds can be considered as a special case of the Grassmannian ones. In particular, for 
m = 2 and n = 1 we obtain the CP^ space where the Berry phases are defined, while for to = n -I- 1 we have the 
identity G{n -\-\,n) « CP". There are two degenerate eigenspaces corresponding to the G{m,n) model that are n 
and m — n dimensional. With the same adiabatic control procedure in the parametric space JVI — G{m, n) non-trivial 
holonomies are generated in both of them, contrary to all previously studied examples. Here, we shall restrict to 
the G(4, 2) case, where the connection components, the corresponding field strengths and a set of holonomies will be 
explicitly given. 
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II. HOLONOMIES AND PHYSICAL MODELS 



A. Berry phases 

The well known Berry phase can be produced by performing loops in the G{2, 1) « CP^ parametric space of 
external control parameters that determine the unitary evolution of two level system as sketched in Figure 1(a). Let 
us present one possible Berry phase implementation in atomic physics. 




(a) (b) 

FIG. 1. Schematic state-structure and interactions parameterized by (a) the G{2,1) and for (b) the G(4, 2) Grassmannian 
manifolds. The arrows represent U{2) transformations between the connected states implemented e.g., by Raman transitions, 
while due to the coset structure there are no transformations between the degenerate states of (b). 



Consider the case of an atom with two non-degenerate ground levels |-|-) and |— ), with corresponding eigenvalues 
= uj/2 and E- ~ —lu/2, and an excited one |e). By performing a Raman adiabatic transfer [10] between levels 
+), |— ) with the help of two detuned laser fields with Rabi frequencies i7+ and fi- and common detuning A it is 
possible to adiabatically eliminate level |e) and create in the basis {|+), |— )} the following unitary transformation 

^ ' \ ie""^sm| cos- 

where = 2\ilcs\t, (j) = argjilcff} and ilcs = {4A). In the above evolution we have neglected the A.C. Stark 

shift effect which can be compensated with properly detuned lasers. In accordance the initial free Hamiltonian of the 
|-f) and |— ) states, given by Hq = j(Jz, is transformed as 

Assume that at time t — we have 9 — (j> — and the initial state of the system is given by |\['±(0)) = |±). Let us 
consider the evolution resulting when a closed path C is spanned very slowly in the parametric plane (0, 0) in time 
interval T. Due to cyclicity and adiabaticity no population has moved from one state to the other at the end of the 
loop, that is, the energy of the system has not changed. Still the states are allowed to obtain a geometrical phase 
factor at the end of this evolution, which can be calculated from the time evolution operator of the system given by 
the Schrodinger equation. A dynamical phase does not appear due to the geodesic character of Raman evolutions. 
In particular, we obtain that |*+(r)) = F;^(C)|-|-) and |*-(r)) = r;^(C)|-) that is the states acquire a geometrical 
phase factor that depend on the spanned loop C, and a connection A. The components of the latter are defined, and 
analytically given for this case, by 

= {±\U^—U\±) = and A^ = {±\U^—Vl\±) = ±ismH . 

In terms of the connection A the phase factors T^{C) = e"^* are given by 

F±(C) = exp / A± = exp // d0dcp Fg%{9, 4>) , (1) 
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where Fg^{9,(f>) = deA^ — d^Ag = ±isin26' is the non-zero component of the field strength associated with A and 
S(C) is the surface on the plane {9, (p) bounded by the loop C. The second step in (1) is due to Stokes theorem. For 
any C we obtain the relation </p_|_ = —if- = JJ^^^^ d9d(f) sin 2^ between the two Berry phases. Note that = 0/2 
where O is the solid angle spanned by the circulation of a unit vector with directions given by the (j)) angles. 
The half factor reflects the spin- 1/2 transformation properties of the employed two atomic levels. This physical 
evolution produces a measurable abelian Berry phase. Their generalization to unitary matrices (holonomies) with the 
employment of degenerate structures is given in the following for the case of the Grassmannian manifold G(4, 2). 

B. Holonomies for the Grassmannian manifold G(4, 2) 

The Grassmannian manifold G(4, 2) corresponds to the parametric space of the Hamiltonian H = UHqU^ where 

Ho = a;/2 diag(l, 1, — 1, — 1) and Li arc U{A) unitary transformations that act non-trivially on the Hamiltonian Hq. 
Clearly, Hq has the following eigenvectors = \l) = (1,0,0,0), \+2) = |2) = (0,1,0,0), = |3) = (0,0,1,0) 
and |— 2) = |4) = (0,0,0,1). The first two states span the degenerate eigenspace S+ with eigenvalue = 
while the last two span the eigenspace S- with eigenvalue = —lo/2. The schematic representation of this model 
is given in Figure 1(b), where the eligible transformations are depicted by arrows. Each arrow correspond to a U{2) 
transformation U{zij) parameterized by a complex number Zij ~ 9ijexpi(j)ij for i = 1,2 and j = 3,4. Hence, the 
real decomposition of the parametric space G(4, 2) is given hy A4 = {9ij,4'ij} = Wij}. Here we adopt the following 
ordering for the general unitary transformation, U{a) = U {zi3)U {zi4)U {z23)U {Z24) ■ 

Now we can define the connection components for each degenerate eigenspace. They are given from the following 
equation 

{AtU^{a\UHa)-^Vliam , (2) 

where the basis vector \a) and belong in the same degenerate eigenspace of Hq- From (2) we see that the matrix 
A^ is anti-hermitian. For {a, 13} = {1, 2} we obtain the elements of the 2x2 matrix A^ , while for {a, (3} = {3, 4} we 
obtain the elements of A~ . The holonomic unitary operator, generated when a closed path, C, is spanned adiabatically 
in the space of the control parameters M., is defined by 

r±(C) = Pexp j> A± . 

The path ordering symbol does not allow to calculate the path integral first and then to exponentiate the result 
as the different components of the connection A do not commute in general. For those cases the simple form of 
Stokes theorem does not apply [11]. Still it is possible to analytically calculate the holonomies which result from 
the commuting components of the connection A in the same way as we did for the abelian case. Alternatively, the 
non-abelian version of Stokes theorem can be employed as in [6,12]. 

Before moving to the analytic calculation of the holonomies, let us see how the transformations U{a) can be 
physically realized for example by an atomic system. According to Figure 1(b) we can generalize the model of the two 
level atom with four levels, pair- wise degenerate and apply laser beams connecting these states via additional exciting 
states by Raman adiabatic transfers. As in the case of the two- level system presented in the previous subsection, a 
U{2) unitary transformation results from each Raman transition. Successive applications of those unitaries are able 
to construct the general unitary transformation U{a) parameterized by the G(4, 2) manifold. A detailed study of 
the generation of the Berry phases with M. = CP^ control manifold with an atomic system manipulated by Raman 
transitions is given in [13]. 

III. CONNECTION AND FIELD STRENGTH COMPONENTS 

A. Connection A 

In this subsection we shall present the connection components related to the parametric space G(4, 2) by employing 
definition (2). In particular, they are 2x2 matrices paired in the following with respect to the S+ and S- degenerate 
eigenspaces they act on. Thus, we obtain 
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/ _gi(013-023) 

At = sin 623 cos cos ^24 I ^ ^ „ . . / , , , , n . 

^4 1^ e-'^-^i^*-^^:*) 2zsin((/)i3-(/.i4 + (/)24- 023) tan 6*14 810 6*24 

(0 _gi(014-013) 

e-'(^"-^") 2z sin(0i4 + <^23 - ./-la - M tan ^23 sin 024 



Af, = sin 

C13 



/ -eii<t>14-4>24) \ _ / \ 

Al^ = sin 024 j , A,^^ " (, j ' 

^■-(ll)^ -..-ill 



The matrix elements of A^^^ are given by 
(^^13)11 = -i sin^ 6*13 cos^ 6*14 , 

1 1 

M+ )i2 = --ie'('^i^-'^==') sin 2013 sin 023 cos 6»i4 cos 024 + 7;ie-'('^"-'^=*)sin2 013 sin 2014 sin 6I24 , 

(vlj^^)22 = i sin^ 013 sin^ 023 cos^ 024 - i sin^ 0i3 sin^ 024 sin^ 0i4 

+ -icos{(j)23 - 013 + 014 - 024) sin 2013 sin 014 sin 023 sin 2024 , 

(^^13)11 = i so? 013 COS^ 023 , 

(A- J12 = - iie-'('^=3-024)sin2 5»^3 sin 2023 sin 024 + ^ie-'('^i3-'^") sin 2013 sin 014 COS 023 COS 024 , 
(^^13) 21 = (^013)12 ' 

(A^^^ )22=i siv? 013 sin^ 023 sin^ 024 — i sin^ 0i3 sin^ 0i4 cos^ 024 

-^icos(-0i4 + 013 + 024 - 023) sin 2013 sin 014 sin 023 sin 2024 , 

.+ ■ 2a ( -1 -e'('^"-'^2*'ctan0i4sin024 \ 

- f^i4 _e-i(<^i4-<^24)ctan0i4sin024 sin^ 024 ) ' 

■ ■ 2n 2n ( ^ ^\ 

^014 = ^ Sl"^ COS 0*24 ( Q ^ I , 

"^^23 = ^ ^^"^^ ^23 COS^ ^24 ^ Q ^ ' 

_. . 2fl 1 e'("^2*-"^23)ctan023sin024 

^^,3 - « Sm 023 ^_i(^,,_^,3)^^^^ _ g.^2 

^^24 = ^si^' ^24 Q ) , A-^^^ = isin^ ^24 J 5 ) 
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B. Field strength F 



In this subsection we shall calculate the field strength F associated with the previous connections A. Its components 
arc given by i^^^ = dfj_Ai, — d^Af^ + [^^, Aj^]. Here, wc shall restrict on the field strength components for which the 
commutator part [A^, A^] is zero for all values of the parameters 9ij and (pij. For them the computation of their 
holonomies is straightforward as we shall see in the following subsection. Let us present first the {6, (j)) field strength 
components. The ones that are obtained from commuting connection components acting on the 5+ and S- eigenspaces 
are given by 

iPL<t>J^^ = ^e^^-^i^-^^^) sin 2013 sin 023 cos 014 sin 024 + ^e-^('^i*-'^^*) sin 2014 sin^ 013 cos 024 , 
(-^^4013)21 = (-^6^4013)12 ' 

(-^^4013 )22 = ~^ sin 013 sin 2024 (sin^ ^14 + sin^ 023) + ^ cos(023 - (^13 + (l>u - ^24) sin 20i3 sin 0i4 sin 023 cos 2024 , 

(^^";40.3)"=O' 

{Fg^4,Ji2 = -le*('^=3-024) sin2 013 sin 2023 cos 024 - ^e'^^^-^i*) sin20i3 sin0i4 cos 023 sin 024 , 

(-^624013)21 = (^624013)12 ' 

{Fe2i<t>i3 )22 = * sin^ 0i3 sin 2024 (sin^ ^23 + sin^ 0i4) - i cos(0i3 - ^14 + (l)2i - ^23) sin 20i3 sin 0i4 sin 023 cos 2024 , 

+ _ f -iie*^^"-'^^^) sin 2014 cos 02 

^ -iie-*('^"-'^2'') sin 2014 cos 024 « sin^ 0i4 sin 2024 

-^^0^4014 = -isin^ 014 sin 2024 q J 

K.<^.. =jsin2 023 sin 2024 ( Q \ 



1 ie'('?*24-</.23) sin 2023 cos 024 



^24023 I iie*(-'^2*+'^23)sin2023COS024 -i sin^ 023 sin 202 



PL<t.. ^ sin 2024 ( n ? ) , F,. = z sin 2024 ° ° 



'24024 \Q I 924024 1 Q ^ 

From the connection components where only the ones that act on S+ commute we obtain the following field strength 
components 

^0^3 013 = 2 ^^'^ ^^^^ 

/ -e*('^i^-'^^3)cos0i4 \ 

^ \^ -e-'C-^ia-fe) cos 014 2 tan 013 sin 023 + 2cos(023 - (i>\3 + (j^u - 024 ) sin 0i4 sin 024 / ' ^ ^ 
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Ft..A..^ = — i cos^ ^24 sin 2^23 
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while -F^3^^^ = -^6^3 ~ ^- From the connection components where only the ones that act on S- commute we obtain 
the following field strength components 

^iii<t'i3 ^ ~ 2 ■^^^^ 

/ -e-'("^i3-"^")cos^23 \ 

^ \^ _e'(<^i3-«ii4) COS6I23 2 tan sin 6114 cos 6124 + 2cos(?!>i3 - (j)i4 + (f>24 - ^23) sin^as sin6i24 / ' ^ ' 

^9T4<Ai4 = ^ cos^ ^24 sin 26114 (01)' 
while -FJ^^023 = Peli4>24 ~ ^' Fin^illy, the {9, 9) field strength components are given by 

^eT,fl,. = sin 6*23 cos 6*14 /. ^ , ^ ^ N.. 



" V e-*('^"-023) 3^6*24 -2zsin((/)i3 - 0i4 + «i'24 - hi) tan 6*14 cos 2^24 7 ' 

p- _ . « n ( -e'('^"-''^i3)sin6i24 

^^13^24 sm cos (^23 g-i(0,4-.^i3) sin 62^ -2i sin(0i4 + <^23 - -^13 - M tan ^23 cos 2^24 ) ' 

/ —Qi(,4>14,-4>24,) \ _ 

^(914024 = - COS 6*24 g_i(0,^_024) J ' ^^^"^13 = 

_ / _gj(024-023) \ 

^023^24 = ^023524 = " COS 6»24 ( q 1 , 

/ _gi(013-<#>23) 

= -COS ^23 COS ^14 COS 6*24 o- • /-J, J, , J, J, \4- a ■ a 

"13023 H<?i3 <P23j 2ism{^i3 — 014 + 024 — fe)tan6'i4sm6'24 

(0 _g«(<^14-<^13) 
e-'(^--^-) 2z sin(0i4 + 023 - <Ai3 - <^24) tan ^23 sin ^24 



C. Holonomies T 



In general, the explicit calculation of the Holonomics of matrix connections A is non-straightforward as they involve 
the path ordering procedure when exponentiating their loop integral. On the other hand it is possible to restrict our 
cyclic evolutions to specific planes (<t, a') that correspond to commuting components A„ and A„i. For those loops we 
can employ Stokes theorem for the abelian theories and obtain 

r±(C) = P exp / yl± = exp (f A^=exp ff dada' F^^, {a, a') , (5) 
Jc Jc J Jt.{C) 

where the rest of the variables are considered constant. The path ordering symbol has been taken out at the second 
step as the connection components on the (a, a') commute with each other. Hence, Stokes theorem can be applied 
straightforwardly as presented in the previous section for the abelian Berry phase. An analytic expression for the 
holonomies can be obtained by exponentiating the 2x2 matrices resulting from the surface integral of the field 
strength components presented in the previous subsection. 
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IV. DISCUSSION 



A theoretical model has been presented for the construction of non-abelian holonomies for the G(4, 2) Grassmannian 
manifold. This is the generalization of the usual abelian Berry phase to the case of a quantum system consisting of a 
doubly degenerate energy eigenspace. The evolution of both degenerate spaces has been presented that are produced 
from the same cyclic adiabatic evolution. The main difference with the abelian case is that it is possible to have 
manipulations of state population in each degenerate space rather than just generation of overall phase factors. This 
can be achieved for example by spanning loops C on the (^24, 4>iz) plane where population can be interchanged in a 
well defined way between the states |1) and |2) as well as between |3) and |4). In contrast, loops on the (^24,024) 
plane contribute only Berry-like phases on the conjugate states |2) and [4) (sec Figure 1). It is worth noticing that 
the holonomic evolution is not producing any correlations between the two different eigenspaces S+ and S- due to 
the cyclicity of the adiabatic procedure. Even though it is possible to have each degenerate eigenspace evolving with 
a different holonomy there is a correspondence between the evolutions as can be easily seen in the previous section 
facilitating eventually their detection in a physical system and the verification of the above results. Indeed, the 
components and F~ of the field strength have a similar functional dependence on the variables {Oij , (j)ij ) and their 
surface integral in (5) bears a common dependence in the area S(C). These holonomic characteristics can be verified 
in the laboratory with present technology by employing an atomic cloud and manipulating the atomic states with 
external laser beams. 
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